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ABSTRACT 

Recent Hubble Space Telescope (HST) observations have revealed that a majority of active galactic 
nuclei (AGN) at z ~ 1 — 3 are resident in isolated disk galaxies, contrary to the usual expectation that 
AGN are triggered by mergers. Here we develop a new test of the cosmic evolution of supermassive 
black holes (SMBHs) in disk galaxies by considering the local population of SMBHs. We show that 
substantial SMBH growth in spiral galaxies is required as disks assemble. SMBHs exhibit a tight 
relation between their mass and the velocity dispersion of the spheroid within which they reside, the 
M. — a e relation. In disk galaxies the bulge is the spheroid of interest. We explore the evolution of 
the M, — cr e relation when bulges form together with SMBHs on the M. — a e relation and then slowly 
reform a disk around them. The formation of the disk compresses the bulge raising its <j e . We present 
evidence for such compression in the form of larger velocity dispersion of classical bulges compared 
with elliptical galaxies at the same mass. This compression leads to an offset in the M, — a e relation 
if it is not accompanied by an increased M. . We quantify the expected offset based on photometric 
data and show that, on average, SMBHs must grow by ~ 50 — 65% just to remain on the M. — a e 
relation. We find no significant offset in the M. — <j e relations of classical bulges and of ellipticals, 
implying that SMBHs have been growing along with disks. Our simulations demonstrate that SMBH 
growth is necessary for the local population of disk galaxies to have remained on the M. — cr e relation. 
Subject headings: black hole physics — galaxies: bulges — galaxies: evolution — galaxies: nuclei — 
cosmology: theory — methods: numerical 



1. INTRODUCTION 

The energetics and demographics of active galactic 
nuclei (AGN ) whi ch are foun d alrea d y by z = 6 
(jZheng et all l2000t iFan et al.l 120001 l200i l200fih . 
be explained by the presence of accre ti ng 
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can be explained by tne presence 01 accre ting su- 
perm as sive black holes ( Lvnden-Belll 119691: 
1981 IChokshi fe Turner! 19921: i Salucci et al.l 
Merritt fc Ferraresd l2001al) . The Hubble Space 



Telescope (HST) has revealed such supermassive 
black holes (SMBHs), with masses in the range 
10 6 — 10 9 M^, in a number of n earby quiescent galaxies 
(|Kormendv fc R ichstonel 119951 ). However the for- 
mation and growth of SMBHs remains something 
of a mystery with a variety of models proposed 
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SMBHs exhibit a number of scaling relations, the 
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tightest of which is the M. — a e relation between their 
mass, M # , and the velocity dispersion, cr e , of the 
spheroids within which they reside. A scaling-relation 
of the f orm logM, = a + /3 1og( g e /200 kms" 1 ) was 
found by lGebhardt et all (|200C0 and lFerrarese fc Merritt] 
(|2000D . Early measurements of the slope /3 varied 
from 4.02 ± 0.32 dTremaine et al.l I2002T) to 4.72 ± 0.36 
(jMerritt fc Ferraresd I2001bf ). More recent measure- 
ments still find a large range of (3 spanning j3 — 
4.24 ± 0.41 (IGultekin et all l2009h to = 5.5 7 ± 0.33 
(jMcConnell et al.ll2011t iMcConnell fc Mall2012l ). Other 
suggested correlations between SMBHs and their host 
bulges include the M, — Lbui or M, — Mbui relations with 



the bulge luminosity o r mass ( Kormendv 
1991 IMagorrian et all [l99l 
Haring fc Rhd 120041 ), and the M. 
bulge Sersic index I 



Richstone 
Marconi fc Hunt! I2003I 



n relatio n with the 

Graham fc Driverl 120071 ) . A three- 
parameter f undamental plane for SMBHs has also been 



suggested (IMarconi fc Huntl 120031; Ide Francesco et al.l 
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200(1 lAlier fc Richstond 120071; iHopk ins et al. I I2007bllal : 
Barwav fc Kembhavill2007[ ). iGrahaml (|2008l ) argued that 
the fundame n tal p lane is caused by barred galaxies. 
iBeifiori et al.l (|20L2T ) found that the fundamental plane 
is strongly dominated by cr e . A relation between M. 
and Mp- „i, the mass of the host galaxy has been sug- 
geste d (jFerraresd 120021 ; IBaes et all 120031 ; IPizzella et at] 
120051 ). Early work found that the same relation is satis- 
fied also by nuclear star clusters (iFerrarese et al.1 120061 : 
iWehner fc Harris! 120061: iRossa et al.l I2006D . which may 
provide a uni fied picture of the growth of central mas- 
sive objects (jMcLaughlin et all 120061 : lHartmann et al.l 
2011). However more recent work has found that nu- 
clear star cl usters and SMBHs follow different scal- 
ing relations (|Erwin fc Gadott i 2012a; EHgheFaDIloil 
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iScott fe "Gr aham 20121). This may possibly be a result of 
the scaling relations being different in late- ty pe galaxies 
(|Greene et al.ll2uTot lErwin k Gadottill2012bft . 

The clues to SMBH growth and formation implied 
by these scaling relations are non-trivial to decipher: 
The sphere of influence of a typical SMBH has a ra- 
dius of a few parsecs, which is some 2-3 orders of mag- 
nitude smaller than the effective radius of a typical 
bulge. What mechanism then gives rise to these scal- 
ing relations? Do SMBHs regulate bulge growth or is 
the growth of SMBHs restricted by the bulge in which 
they reside? Examples of the latter are scenarios in 
which gas accr etion onto the SMBH is regulated by 
star fo rmation (|Burkert k Silkl 120011 : iKazantzidis et all 
[2001 iZhengetal] |2009ft or by stellar feeding of 
SMB H accretion disks ([Miralda-Escude k Kollmeierl 
[2001 . AGN feedback via heating, pressure-driven 
winds or ionization typically gives rise to scenarios 
in which SMBHS r egulate their own or b ulge growth 
(ISilk k Reesl Jl998) IWvithe k Loebl [200l lKrn¥l200l 



Murray et alfeoOsflDi Matteo et al.H2005HSazonov et alj 



20051: lYounger et all 12008ft . High velocity outflows that 



may be associated with such AGN fee dback have been 
obser ved in Seyfert 1 galaxies (e.g. iCrenshaw et al.l 
11999ft . In addition, in semi-analytic models, AGN 
feedback is often invoked to exp lain the high mass 
end of the lumin osity function (|Croton et al.l 120061 : 
iBower et al.l [2006). Alternatively, collapse mod- 
els in which the M. — a e relation is an indirect 
consequenc e of unrelated proces s es have also been 
proposed ( Hachuc rt k Kauffmannl 120001: lAdams et al] 



20011 1200a iPend l2007t iVolonteri k Nataraianl 120091: 
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One of the characteristics of the scaling relations 
is that M, correlates with the properties of the host 
spheroid. In disk galaxies this is the bulge component. 
Bulges in disk galaxies come in two types: "classical" 
and "pseudo" bulges, with mixed typ es als o possible (e.g. 
Erwin et al] 12003: jDebattista et all 120051: lAthanassoulal 
20051 : iNowak et al.1 120101 )! Classical bulges are believed 
to form via merging of sub-gal a ctic clumps, satellites 
and clusters (lEggen et al.l 119621: iTremaine et al.l 119751: 



Searle k Zin^[l978tlKauffmann et al.lll99aiBaugb. et al] 
19961 : Ivan den Boschlll998ir In essence classical bulges 



are ell iptical galaxies around which a disk has reformed 
(e.g. iSteinmetz k Navarro! 120021 ) although contin- 
ued late growth of clas sical bulges is also possible 
(e.g. iHopkins et al.l [20To| ). Pseudo bulges instead are 
formed by the secular evo lution of disk structure , 
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reviewed the observational evidence for pseudo bulge 
formation. In contrast to pseudo bulges, classical 
bulges form early, predating the formation of the 
disk. The difference between classical and pseudo 
bulges is reflected als o in th eir SMBH demographics. 
lGadotti"fc Kauffm annl (|2009[ ) estimate that classical 
bulges account for 41% of the black hole mass in 
the local unive r se, wh ile pseudo bulges host only 4%. 
IGultekin et al.l ([2009ft find that the classical bulges 
(including elliptical galaxies) follow the same M. — a e 



relation as the general populatio n, with a scatter of 
0.45 ± 0.066. Instead [Hyf (120081 ) found that pseudo 
bulges have an M, — o~ e relation with the s ame slope as, 
but lo wer zero-point than, classical bulges. iGreene et al.1 
(|2010ft showed that the SMBHs of late-type galaxies, 
which predominantly contain pseudo bulges, scatter 
below the M. — a e relation. 

A number of studies using simulations have explored 
the evolution of the M, — o> relation during hierarchi- 
cal mergi ng (e.g. IKazantzidis et al.l 120051: iSpringel et all 
2005alrbT lYounger et alj 120081 : Uohansson et al.l 120091 : 
Robertson et al.l 12006ft . However, recent HST observa- 



tions have shown that a surprisingly substantial frac- 
tion of AGN activity at high redshifts is associated 
with isolated dis k galax ies, rather than with mergers. 
ISchawinski et al.l ([20111 ) show that ~ 80% of X-ray- 
selected AGN at z = 1.5 — 3 are in low Sersic-index galax- 
ies, indicative of disks. They find that moderate luminos- 
ity AGN hosts at z ~ 2 are similar to those at z ~ 0. Ex- 
cluding the high luminosity qua sars, which are triggered 
by mergers (|Treister et all2"010ft , they estimate that 23 — 
40% of SMBH growth occurs in intermediate brightness 
Seyfert AGN. The X-ray-selected samp le of moderate- 
lumin osity AGN at 1.5 < z < 2.5 of iCisternas et afl 
(|2011ft consists of more than 50% disk galaxies, with on- 
going mergers e vident no more f r equen tly than in non- 
active galaxies. ISchawinski et al.l ([2012ft show that even 
heavily obscured quasars are hosted largely by disks, not 
by mergers. Studies of star-formation using Herschel 
find that the specific star formation rates of X-ray se- 
lected AGN hosts are no different from those of inactive 
galaxies, also indicating that AGN host s are not under- 
going fundamentally different behaviors ([Mullanev et al] 
l2012bl lal) . Using m ultiwavelength survey s of AGN across 
redshifts < z < 3 lTreister etaTI (|2012ft found that only 
the most luminous AGN phases are connected to major 
mergers, the rest being driven by secular processes. The 
merger driven AGN activity accounts for only ~ 10% of 
AGN. The "anti-hierarchical" nature of galaxy and AGN 
growth — both the largest galaxies ([Bower et alj 119921 : 
iThomas et~aLl 12005b [Nelan et al] 12005ft and the bright- 
est AGN ([Ueda et al.ll2003l : iHasinger et al .112005ft form at 
high redshift whereas lower mass galaxies and moderate 
luminosity AGN peak at lower redshifts — also hints that 
internal evolution rather than mergers is the main driver 
of SMBH growth. AGN activity continues to be domi- 
nated by disk galaxies down to the present: since z ~ 1 
more than 85% of AGN activity is hosted in galaxies wit h 
no evidence of recent mergers ([Kocevski et al] 12012ft . 
Lastly, the presence of AGN in bulgclcss galaxies, which 
are thought to not have experienced much hierarchical 
merging, provides further evidence that internal evolu- 
tion is capable of driving SMB H growth ([Simmons et alJ 
l2T^lArava Salvo et al.ll2012[) . 

This paper introduces a novel approach to exploring 
the origin of the M, — <r e relation. We study the con- 
sequences of disk regrowth for the M. — a e relation of 
classical bulges under the assumption that a classical 
bulge forms with a SMBH satisfying the M, — a e rela- 
tion and then a disk reassembles aro und it. Growth o f 
the disk then compresses the bulge (jAnd rcdakis 1998ft . 
Since o~ e is not an adiabatic invariant, compression leads 
to its evolution, which we quantify here. We study the ef- 
fect of this evolution on the M, — <j e relation. If SMBHs 
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remain on the M, — a e relation then this implies that 
SMBH growth is governed by the potential (as charac- 
terized by cr e ) within which they sit, which is most likely 
if AGN feedback regulates SMBH growth. If instead we 
find that bulges evolve away from the M. — a e relation, 
with the SMBHs retaining a memory of the bulge within 
which they formed, then this implies that bulge growth 
is limited by the SMBH, as would happen if AGN feed- 
back quenches star formation in the bulge. The paper 
is organized as follows. Section [2] describes the simula- 
tion methods used in this paper. Section [3] presents the 
evolution of a e caused by disk (re-)assembly and derives 
a photometric estimate for the increase in er e . In Sec- 
tion [4] we predict the consequences of bulge compression 
for the M, — a e relation. We find that the main effect 
is a shift to lower mass in the zero-point of the relation. 
Then in Section[5]wc test this prediction on observational 
data. We find no evidence for such a shift, indicating that 
SMBHs have grown along with disks. We show that the 
degree by which SMBHs must have grown is consistent 
with the new HST estimates. Section [6] sums up our 
results. 

2. NUMERICAL SIMULATIONS 

We construct initially spherically symmetric two- 
component galaxy models consisting of a stellar bulge 
embedded in an extended dark matter (DM) halo. 
For the DM compone nt, we consi d er th e cuspy, 
cosmologicallv-motivated lNavarro et al.l (|1996l hereafter 
NFW) density profile given by 



pr>M(r) 



(r/r s ){l + r/r s y 



(r < ryir), (1) 



where p s is a characteristic inner density, r s denotes 
the scale radius of the density profile defined as the 
distance from the center where the logarithmic slope, 
dlnp(r)/dlnr, is equal to —2, and r v i r is the virial ra- 
dius defined as the radius enclosing an average density 
equal to the virial overdensity times the critical density 
for a flat universe. We adopt the ACDM concordance 
cosmology with Q m = 0.3, = 0.7, and h = 0.7, and 
assume z = 0. T he virial overdensity is then equal to 
A vir ~ 103.5 fe.g- ILacev fc Coldll9m 

The NFW density profile is formally infinite in extent 
with a cumulative mass that diverges as r — > oo. In order 
to keep the total mass finite, we implement an exponen- 
tial cutoff which sets in at the virial radius and turns 
off the profile on a scale r^ ecsiy . The truncation scale 
''decay is a free parameter and controls the sharpness of 
the transition. Explicitly, we model the density profile 
beyond r v ; r by 



p(r) 



P s 



c(l + c) 2 





r - r vir " 


exp 




\ ' vir / 


^decay 



(r > r V ir), 



(2) 

where c = r v i T /r s is the concentration parameter and k 
is fixed by the requirement that dlnp(r)/dlnr is con- 
tinuous at r vir . This procedure is necessary because 
sharp tr uncations result in mode ls that are not in equi- 
librium fKazantzid is et al.l |2004j) . For the purposes of 
the present study, we adopt a concentration parameter 
c = 10, app ropriate for Milky W ay Galaxy-sized dark 
matter halos (|Bullock et al.|[200ll) , and a truncation scale 

'"decay — 0.l7*g. 



For the spatial distribution of the b ulge compon ent, 
we adopt the de-projected Sersic law (jSersid |1968[ ) of 
iSimonneau fc Pradal {2001 : 



P(s) = Po 



where 



exp 



-ks" (1 — x 2 ) 



Po 



l-(l-x 2 ) — 
k Y. 2 1 



da: 



7r R e o n - 1 



(n>l), 
(3) 

(4) 



In the above equations, n denotes the Sersic index, R e Q 
is the effective radius, i.e. the radius that encloses half 
the total projected luminosity, T,q is the central value of 
the projected mass profile, and s = r/i? e ,o- For n > 1, 
k can be estimated (with an e rror smaller t han 0.1%) by 
the relation k = 2n - 0.324 (jCiottil fl99l[) . Our initial 
b ulge has Sersic index n = 4, i.e. it is characterized by 
a Ide Vaucouleursl (|1948f ) profile. For the specific galaxy 



model we consider, the ratio between the mass of the 
bulge and the virial mass of the halo is equal to 8 x 10~ 3 , 
while the ratio between the bulge effective radius and the 
halo scale radius, R e fi/r s = 0.02. Our standard scaling 
has i? e ,o = 500 pc and bulge mass Mb = 8 x 10 9 M , 
leading to r vir = 270 kpc and M vir = 10 12 M Q . Because 
we do not consider non-gravitational processes such as 
gaseous dissipation, the scale-free nature of gravity al- 
lows the rescaling of our models. 

Monte-Carlo realizations of the iV-body galaxy model 
are constructed accordin g to the procedure described in 
iKazantzidis et ail (|2004f) , which is based on sampling the 
exact phase-space distribution function (DF). Under the 
assumption of isotropy, the DF of each component de- 
pends only on the binding energy per unit mass E: 



ME) 



1 



787T 2 



d* 



1 



d^ 2 y/E- * y/E 



dpi 
d* 



(5) 

where pi is the density profile of component i and ^(r) = 
ipDM(r) + '/'stars (?") is the total relative gravitational po- 
tential. Note that the second term on the right-hand side 
in Eqn. [5] vanishes for any sensible behavior of ^(r) and 
Pi{r) as r —t oo. 

The system generated this way needs to be softened; 
the gravitational softening lengths are set to e = 15 pc 
for all particles (including the disk particles described be- 
low) in all runs. Since softening the pote ntial is equiva - 
lent to smoothing the density distribution (|Barnesll2012[ ) , 
the initial conditions set up without softening are not 
a perfect equilibrium. We therefore relax the initial 
bulge+halo system for 250 Myr before we start growing 
the disk. During this period the bulge settles to a new 
equilibrium. For the remainder of this paper we refer to 
this relaxed model as the initial conditions. 

After the bulge+halo system has reached equilibrium, 
we investigate its response to the growth of various exter- 
nal disk fields. The growing disks follow an exponential 
distribution in cylindrical radius R, and th eir structure 
is modeled as (lSpitzerlll94l iFreemanll 19701 ): 



Pd(R,z,t) 



m d (t) 
8irz d R 2 , 



exp 
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Rd, 



4 



Dcbattista et al. 



where m</, Rd-, and Zd denote the mass, radial scale- 
length, and vertical scale-height of the disk, respec- 
tively. Except for one model, we use Zd = 0.15 Rd in 
all experiments, a choice which is consistent with ob- 
servations of external galaxies (Ivan der Kruit fc Sear lei 
Il98l Ide Grilsfc 7 van der KruTtlll996l) .~ However, the ob- 
served scatter in Zd/Rd is quite substantial, reaching 
to values as small as 0.05. We show below that thin- 
ner disks lead to even stronger compression. Thus our 
assumption of Zd = 0.15 Rd is conservative. We im- 
plicitly assume that the cla ssical bulge is fully f ormed 
at the last major merger. IHernquist fc Mihoa (|1995D 
showed that minor mergers drive gas to small radii lead- 
ing to a burst of star formation and bulge growth. How- 
ever dissipation is now thought to largely give rise to 
pseudo, not classical, bulges. The origi n of bulges in 
high m ass galaxies remains contentious. IWeinzirl et al.1 
(|2009f) used bulge+disk+bar decompositions to argue 
that mergers cannot account f or the majority of bu lges 
in current high mass galaxies. iHopkins et al.l (|2010() in- 
stead argued that major mergers dominate the forma- 
tion of these bulges, with minor mergers contributing 
another ~ 30%. Nonetheless a separation between clas- 
sical and pseudo bulges seems to be well established, 
with observational evidence indicating that properties 
such as morphologies, star formation rates and cor- 
relations with d isk properties, includi n g color, change 
across n ~ 2 (IDrorv fe Fisher! 120071: [Fisher fc Drorvl 
l200aiFisher et al.H2009UFisher fe Drorvll2010D . 

Each growing-disk simulation is performed by grow- 
ing linearly over time the mass of an initially mass- 
less Monte Carlo particle realization of the desired disk 
model: rrid{t) = (t/r)Md. In all simulations, we set 
M d = 1-6 x 1O U M and r = 2 Gyr. During the exper- 
iments the disks are held rigid with their particles fixed 
in place, while both the bulge and halo particles are live, 
allowing them to remain in equilibrium as the disk mass 
grows. Throughout the experiments, all other properties 
of the growing disks (e.g., scale- lengths, scale- heights) 
are kept const ant. Additional details of this technique 
can be fou nd inlDebattista et all (120081) , iVillalobos et al.l 
(|2010D and lKazantzidis et al.l (|2010l) . Our simulations do 
not include a SMBH since the sphere of influence of the 
SMBH that would correspond to the velocity dispersion 
of the initial bulge, = GM./er e 2 is only ~ 4 pc, which 
is considerably smaller than i? c .o- 

The initial conditions of the two-component galaxy 
contain a total of 4.4 million particles (4 x 10 6 dark 
matter particles and 4 x 10 5 bulge particles). The disk 
is modeled with a further 4 x 10 5 particles . Parti- 
cles a re set up using a quiet start procedure (jSellwoodl 
119831) that ensures that all components have zero net 
momentum. We set up particles in groups of four: 
the first particle has (x,y, z,v x ,v y ,v z ) while the rest 
have (-x,-y,z,-v x ,-v y ,v z ), (-y,x,-z,-v y ,v m ,-v z ) 
and (y, -x, -z, v y , -v x , -v z ). 

We run 5 simulations, with varying ratio of disk scale- 
length to initial bulge effective-radius: R,x/Re,o = 0.5, 1, 
2, 5, and 10. We save outputs at 25 Myr intervals corre- 
sponding to increments SMd = 0.25M;,. Table [T] provides 
a summary of a representative subset of the outputs. 

All numerical simulations a re carried ou t with the par- 
allel A^-body code PKDGRAV (jStadelll200lT ). In all exper- 



iments, we set the base timestep At = 1.25 Myr with 
timesteps refined such that St = At/2 P < ^(e/a) 1 / 2 , 
where e is the softening and a is the acceleration at a 
particle's current position, with rung number p as large 
as 29 allowed. For all simulations we set r\ = 0.02 and 
use an opening angle of the treecode 8 = 0.7. In the 
Rd/Rc,a = 1 model, timesteps for particles get as small 
as 2 -10 of the base timestep, (i.e. 1220 years). 
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0.5 


10 


10.04 


113.4 ±0.9 


0.9 


1 


10 


10.14 


113. 7± 0.5 


0.9 


5 


10 


10.68 


116.9 ±0.8 


1.1 


10 


10 


11.22 


119.4 ±0.8 


1.2 


15 


10 


11.56 


122.2 ± 1.3 


1.3 


20 


10 


11.80 


123.5 ± 1.3 


1.3 



TABLE 1 

A REPRESENTATIVE SAMPLING; OF THE RESULTS AT DIFFERENT 
TIMES (CORRESPONDING TO DIFFERENT D/B) OF THE SIMULATIONS 
PRESENTED IN THIS PAPER. THE FIRST TWO COLUMNS LIST THE 
'INPUT' PARAMETERS, WHERE D/B IS THE DISK-TO-BULGE MASS 
RATIO AND Rd/Rc,0 IS THE RATIO OF DISK SCALE-LENGTH TO 
INITIAL BULGE EFFECTIVE RADIUS. THE LAST THREE COLUMNS LIST 

THE OUTPUT PARAMETERS OF THE SIMULATIONS. Rd/Re IS THE 
RATIO OF DISK SCALE-LENGTH TO FINAL BULGE EFFECTIVE RADIUS, 

<T e IS THE VELOCITY DISPERSION OF THE BULGE AND T. IS THE 
RATIO OF FINAL TO INITIAL M. ASSUMING THE SYSTEM STARTS AND 
ENDS ON THE SAME M. — <T e RELATION WITH SLOPE j3 = 4. BOTH 
Rd/Rc AND <J e INCLUDE THE EFFECTS OF RELAXATION OF THE 
INITIAL CONDITIONS. THE FIRST ROW CORRESPONDS TO THE 
INITIAL CONDITIONS BEFORE THE DISK IS GROWN. 



3. EVOLUTION OF VELOCITY DISPERSION 

Growth of the disk compresses the bulge and raises 
its velocity dispersion everywhere. In Figure [1] we plot 
examples of this evolution. From an initial value of 
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~ 100 kms 1 , the bulge in the i?d/-R .o = 1 case attains 
values of oy ~ 300 kms -1 . At each output we measure 
the bulge effective radius, R e , by computing the circu- 
lar projected radius containing 50% of the bulge mass at 
inclinations i = 0°,60° and 90°. A detailed analysis of 
the evolution of the structural parameters will be pre- 
sented elsewhere. We then measured a e in slits as the 
root-mean-square (rms) of the line-of-sight velocity of all 
particles within R e 

R 

2 Eiefl. ratios Jo c I(R)(vf os + (T? os )dR — 

a e = ^ ' = n , {<) 

^ieR, m i ff* I{R)dR 

where vi os is the mean line-of-sight velocity and ai os is 
the line-of-sight velocity dispersion along the slit, which 
is placed along the major (i.e. inclination) axis. We re- 
peat the measurements of a e for the same set of incli- 
nations (i = 0°, 60° and 90°) and use the average of 
these measurements for er e . Because <r e is not identical 
for all viewing orientations, we use the largest difference 
between the average o~ e and the individual values as an 
estimate of the uncertainty on o~ e . We use slits of width 
20 pc, which is ~ 4% of R e o and ~ 16% of the smallest 
R e attained by the simulations, ~ 0.14 kpc. The change 
in <j e is not a result of different radial sampling caused by 
the change in R e , but is genuinely caused by evolution of 
the velocity dispersion, as can be seen in Figure [TJ We 
also measure as, the velocity dispersion in apertures of 
Re/8, similarly by restricting the sum in Eqn. [7] to that 
smaller radius. 




0.5 1 1.5 



r / R e.O 

Fig. 1. — Evolution of ay, the radial component of the velocity 
dispersion in spherical coordinates, for the bulge particles in the 
Rd/Rc,o = 1 case. The profiles are at D/B = 0, 5, 10, 15 and 20, in 
increasing order. The arrows indicate R e in each case. 

The main simulation results are presented in Table [T] 
and plotted in Figure [2] In order t o provide a fitting 
formula to these values, we note that iWolf et al.l (|2010h 
find Mi/2 — 4G~ 1 i? c (<ti os 2 ) for pressure-supported sys- 
tems, where M1/2 is the mass within the half-mass ra- 
dius and (c; os 2 \ is the luminosity-weighted square of the 
line-of-sight velocity dispersion over the entire system. 



We therefore expect a to scale as 

( mosj{R f ) \ 2 = ( MbjjRf) + M d (R f ) \ (Ro\ ( ) 
\o-ios,o( R o) J V M h0 (R ) J \RfJ 1 ' 

where subscripts / and indicate final and initial values 
and Mf,(R) and Md(R) indicate bulge and disk masses 
within radius R. Integrating over radius, we expect 

£_( 1+ 4_( 1+ £).-*/«.]§y 

= T (9) 

where o e $ is the initial dispersion of the bulge. We stress 
that JF is defined to be a photometric, not kinematic, 
quantity. Note that at fixed disk-to-bulge mass ratio, 
D/B, the minimum in T occurs at R c jR^ = 0, i.e. as 
i?d — > 00. In deriving Eqn. [9] from the more general 
Eqn. [S] we have assumed that the disk is exponential; 
this is true for our simulations, but need not b e the case 
in nature (e.g. IBoker et al.ll200l IDuttonl 12009). If 7 = 1 
then the term in the outer brackets on the right hand 
side of Eqn. |H] is merely the ratio of final (bulgc+disk) 
mass to initial (bulge only) mass within R e . (We ig- 
nore the dark matter halo in this calculation since the 
dark-to-bulge mass fraction within R c ^ is less than 2%.) 
We have also assumed that the disk scale-height is small 
compared to the effective radius of the bulge, but the 
factor 7 is introduced to account for some of deviations 
resulting from this assumption. We have neglected the 
compression of the bulge in deriving this expression, i.e. 
we assume that R e = i? .o- which Table Q] clearly shows 
is not the case. We fold the uncertainty resulting from 
this assumption into the free parameter 6, which would 
be 0.5 if R c = R c0 . The best fit for 2 < R d /R c < 9, 
D/B > and 1 <V e /cr e ,o < 2 is 7 = 0.3 and S = 1.76 
which gives a x 2 = 172 for 196 data points. The best fit 
value with 7 = 1 is 8 = 0.64 with x 2 = 562. Likewise, we 
fitted best-fit parameters for ag obtaining 7 = 0.02 and 
S = 15.92 with a x 2 = 719. These best fits are shown in 
Figure [U Because the compression of the bulge depends 
only on the ratios of bulge-to-disk masses and sizes, Eqn. 
[S] remains true for any galaxy or bulge mass, i.e. for any 

c e , ci- 
In order to check whether the bulge and halo respond 
adiabatically to the growth of the disk we slowly evapo- 
rated the disk from the final state of the simulation with 
Rd/Rc.o = 1- The results of this test are shown in Figure 
[3] and indicate that indeed the response is adiabatic to 
good approximation. 

3.1. Other dependencies and sources of scatter 

Eqn. [5] allows us to estimate the increase in o~ e for a 
bulge given a galaxy's photometric decomposition. We 
now explore the amount of scatter that can occur in these 
estimates for a given observed density distribution. Some 
of the effects we consider here will compress the bulge to 
a different extent, so both a e and R c will change; it may 
happen however that Eqn. [5] still provides a good ap- 
proximation to the evolution of a e . Because we use Eqn. 
[S]to estimate the offset of galaxies from the M. — a e rela- 
tion, we are primarily interested in those changes which 
Eqn. [9] does not reproduce, and we consider this to be 
the scatter of interest here. 
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Fig. 2. — Left: <r e as a function of D/B (top) and Rd/Rc,o (bottom). Gray, blue, green, cyan and red represent, respectively, Rd/Rc,o 
= 0.5, 1, 2, 5 and 10 in the top panel and D/B = 1, 5, 10, 15 and 20 in the bottom panel. The initial conditions are indicated by the 
black point in the top panel. Right: <r e (top) and erg (bottom) as functions of R d /R e . Colors are as in the bottom-left panel. The solid 
lines show our fitting function, Eqn. [9] In each panel, the dotted horizontal lines show contours of constant (<t/<to) 4 , with values indicated 
above each contour. 



The (re)-assembly of disk galaxies is not necessarily 
a slow, adiabatic process. A possible source of scatter 
might therefore be due to disks growing more rapidly 
than assumed here. In order to test what the effect of 
faster disk growth may be, we grow the disk ten times 
faster, i.e. within 200 Myr. The results are shown in the 
top panel of Figure HI The effect of the different growth 
rate is negligible. 

At a given mass, the concen tration of dark matte r ha- 
los can vary substantially (e.g. lWechsler et al"1l2002f) . We 
explore what effect this might have on a e by re-running 
our simulation with a halo having c = 20. The bot- 
tom panel of Figure 0] shows that a e is barely changed, 
undoubtedly because the galaxy is baryon-dominatcd in 
the bulge region. Reasonable variations in halo concen- 
tration therefore do not produce any significant scatter 
in the M, — a e relation. 

All simulations above used z& = 0.15i?d- In the top 
panel of Figure [5] we show the effect of halving Zd . This 
increases the final a e to 268.7±13.3kms~ 1 , an increase by 



6.8% while R e decreases by ~ 3.3%. The bottom panel 
of Figure [5] shows that, taken together, these differences 
lower the quality of the fit of Eqn. |H] for Zd = 0.075i?d 
compared with that for the standard Zd = 0.15i?d, al- 
though all cases still have errors of less than 15%. At 
Oe/oe.o = 1.8 the maximum error is about 12%, which 
we adopt as our estimate for the scatter due to disk thick- 
ness. 

Another source of scatter comes from the contamina- 
tion of measured bulge kinematics by the kinematics of 
the disk, which is almost inevitable in real galaxies. Ex- 
ploring this effect requires that we set up equilibrium 
kinematics for the disks in the Rd/R E ,o = 1 model at 
various values of D/B. We set the kinematics of the 
disks to give constant Toomr e-Q = 1.5, as described in 
iDebattista fc S cllwood (2000). For this we calculate the 
potential using a hybrid polar-grid code with the disk on 
a cylindrical gri d and the bulgc+halo on a spherical grid 
(jSellwoodl 120031 ) . Figure |6] shows the effect of disk con- 
tamination: changes in a e can be either positive or neg- 
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Fig. 3. — Evolution of <r e in the Rd/R c ,o = 1 model as the 
disk is grown (red shaded region) and subsequently evaporated 
(black points with error bars) . The dotted horizontal line indicates 
oe/oefi = 1.8 for which Y. = 10.5 if /3 = 4. 
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Fig. 4. — Evolution of cr e in the Rd/R c ,o = 1 model for different 
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black points with error bars show the variant simulation results. 
The dotted horizontal lines indicate a- e /fT e ,o = 1-8 for which T. = 
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Fig. 5. — Comparision of the evolution of <r e in the R^/R c o = 1 
model with a thinner disc Zd = 0.075-Rd instead of the fiducial 
Zd = 0.15i?d. Top: The shaded red region represents the evolution 
in the fiducial case, while the black points with error bars show the 
thinner disk. The dotted horizontal line indicates a e /cr ey o = 1.8 for 
which r. = 10.5. Bottom: A comparison of the fit of Eqn. [9] for 
the fiducial case (black open squares) and with Zd = 0.075i?d ( re d 
filled circles). The dotted vertical line indicates cr e /rr B o = 1.8. 

ative, but generally \a e {B) - cr e {B + D)\/a e (B) < 25%. 
The error increases with a e , which is a result of the in- 
creasing D/B. An independent analysis of the effect of 
disk contamination on a e by Hartmann et al. (in prepa- 
ration) also finds fractional changes < 25%. The open 
stars in Figure [6] show the inclination-averaged values 
of o e (B) versus <r e (B + D); the differences between the 
means are generally less than 20%. However on average 
<j e (B + D) is systematically larger than a e (B). 

Based on these tests we conclude that there is < 30% 
uncertainty in the degree to which classical bulges are 
compressed. 

3.2. Observational evidence for bulge compression 

We now present evidence that bulge compression as- 
sociated with disk rcgrowth has occurred in nature by 
comparing the properties of classical bulges and ellipti- 
cal galaxies. This requires a large sample of disk galax- 
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Fig. 6. — The effect of including disk stars in the measurement 
of <r e . <r e (B) represents the measurement from bulge particles only 
while cr e (B+D) includes disk particles in the measurement. Dotted 
lines have constant slope, as indicated along each line. The different 
points represent the effect of disk contamination as D/B increases 
from 2.5 to 20. Different filled symbols correspond to different 
galaxy inclinations, as indicated. Open (red) stars correspond to 
inclination-averaged values. 



ies with bulge+disk decompositions. iGadottil (|2009() pre- 
sented a detailed structural analysis of nearly 1000 galax- 
ies fro m the Sloan Digital Sky Survey (SDSS) (jYork et al.l 
I2000T ) . classifying them into ellipticals or disks, distin- 
guishing the l atter by whether they host classical or 
pscudo bulges. iGadotti fc Kauffmannl (|2009t) present the 
velocity dispersions within R e /8, hereafter as, for a frac- 
tion of these galaxies. We use these data to compare the 
distributions of as for ellipticals and classical bulges. The 
sample contains 196 elliptical galaxi es and 176 unb arred 
classical bulges with kinematic data. IGadottil (|2009D clas- 
sifies the bulges based on the Kormendy (me an effective 
surfac e brightness (p, e ) versus R c ) relation (jKormendvl 
11977ft . iFisher fc Drorvl ([20081 ) identify Sersic index n = 2 
as the dividing line between pseudo and classical bulges, 
with the latter having n > 2. We apply this additional 
criterion to the sample, which leaves 166 galaxies as our 
final sample of unbarred, classical bulges. 

Figure [7J plots the distribution of galaxies in the 
ag-Mbui plane for ellipticals, observed unbarred clas- 
sical bulges, and the same classical bulges if they 
are decompressed using Eqn. |H1 We obtain Mbuli 
D/B and i? e /-Rd from th e the exponenti al disk+ Sersic 
bulge decompositions of IGadottil (|2009[ ) and <r 8 from 
IGadotti &: Kauffmannl ([2009D . We decompress to obtain 
Us, o using Eqn. [9] fitted to erg. 

The line in the top panels of Figure [7J shows the fit 
to the ellipticals: crg/kms" 1 = (M bu i/3051M Q ) a30783 . 
The observed classical bulges have larger as, on average, 
than the ellipticals at a given Mb u i- When we decom- 
press the bulges their offset from the elliptical relation is 
significantly reduced, as can be seen in the bottom pan- 
els of Figure [7j The distributions of the residuals from 
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Fi g. 7. — Distribu tion s of ellipticals and class i cal b ulges, taken 
from IGadottil j20QgT) and IGadotti & Kauffmannl ^2009) . in the trg- 
Mbui plane. Left: Elliptical galaxies. Center: Classical bulges in 
unbarred galaxies. Right: The same sample as in the central panel, 
but using the decompressed values for ag obtained by applying 
Eqn. [9] The (red) solid lines in each of the upper panels show 
the best fit to the elliptical galaxies. The bottom panels show 
the residuals for each sample from this best fit to the ellipticals: 
A log erg = log erg — log erg fl t , where erg g t is the erg value from the 
fit to the ellipticals of a given mass. 

the fit to ellipticals are shown in Figure [5] The means of 
the residuals are 0.06 dex and 0.04 dex for the observed 
and decompressed classical bulges, respectively. A two- 
sample unbinned K-S test comparing the ellipticals and 
bulges shows that the probability that the residuals are 
drawn from the same distribution is 3 x 10~ 9 for the ob- 
served bulges and a much larger, though still formally 
small, 3 x 10~ 6 for the decompressed bulges. In Figure 
[5] we plot the distribution of the bulges and ellipticals i n 
the as~R e plane. As was also found by IGadottil (|2009t) . 
the observed classical bulges are offset to larger as and 
smaller R c relative to the ellipticals in this projection of 
the Fundamental Plane, as expected if bulges are com- 
pressed by disks. 

While Figures [7J and [S] do provide evidence for bulge 
compression, they also show that the as is even larger 
than predicted by our simple model. Possibly this is 
because disks are more concen trated than expon ential 
at the c enter as proposed by IBoker et all ([2003D and 
IDuttonl (|2009f ). Alternatively, contamination of as by 
the disk, or differences in formation histories, could be 
to blame for (part of) the offset between ellipticals and 
classical bulges. We note that if disks are more con- 
centrated than exponential then their effect is to further 
compress bulges and we have under-estimated the evo- 
lution of a e . Exploration of this issue is deferred to a 
future publication. 

4. CONSEQUENCES FOR THE M. - cr e RELATION 

The steepness of the M. — a e relation implies that the 
maximum factor of ~ 3 increase in a e obtained in the 
simulations would require a factor of ~ 80 increase in 
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Fig. 8. — Distributions of the residuals shown in the lower panels 
of Fig. [7] The solid (black) line is for the ellipticals, the dashed 
(red) line is for observed classical bulges (as described in the text) 
and the dot-dashed (blue) line is for the classical bulges decom- 
pressed using Eqn. [5] 
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Fig. 9. — The ag-R e p r ojecti on of the fundamental plane for 
the Gadotti & Kauffmann (2009) sample. The (black) filled circles 
show elliptical galaxies while the (red) open circles show the ob- 
served unbarred classical bulges. The diagonal (blue) line shows 
the evolution of the model in the simulation with R d /R Ct o — 2, 
assuming ii c ,o = 2 kpc and <7 e ,o = 100 km s . The star symbols 
correspond to the system at D/B = (bottom right), 5, 10, 15 and 
20 (top left). 

M, for the SMBH to remain on the relation, more than 
6 x larger than the factor of 20 by which the stellar mass 
grew. We define the factor by which M. must grow 
to remain on the M. — a e relation, T. = M,j/M,fi = 
(cr e /cT e .o)' 3 , where subscripts and / indicate initial and 
final values. The dotted contours in Figure [2] indicate 
r, assuming j3 = 4. Values of T, for the simulations 
are listed in Table Q] In general T. > 3 requires that 
Rd/Re < 5 and D/B > 2. 
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Fig. 10. — Dependence of ag on the parameters of bulge +disc de- 
compo sitions for unbarred classical bulges in the sample of lGadottH 

4.1. Evolution of slope and zero-point 

In Figure[lU]we plot the dist ribution of cr s as a function 
of D/B and of R d /R c for the IGadottil ([20091) sample; a 
weak correlation between D/B and o~$ is present (Spear- 



man r 



sp 



-0.22, Kendall r = —0.15), which is statisti- 
cally significant at less than 3c. The correlation between 
Rd/R c and erg is even weaker (Spearman r sp — —0.10, 
Kendall r = —0.07). We therefore neglect these weak 
correlations. If cr e /a e .o — J 7 , and the average a e /o~efi at 
a given er e u is (J 7 ), then neglecting these correlations im- 
plies that (J 7 ) is independent of cr ej o- Then wc can write 
the M. — a e relation for compressed bulges, if M. does 
not change while disks grow, as 



logM. =a-p\og (J 7 ) + f3 log a e 



(10) 



i. e. the slope of the relation remains /3 but the zero-point 
changes by 

5a = -p\og{T). (11) 

Because J- > 1 (compression can only increase er e ) Eqn. 
ITTI implies that Sa < 0. Therefore, if M. does not 
grow during disk formation, the M, — a e relation of com- 
pressed bulges will be parallel to, but offset below, the 
M. — a e relation for elliptical galaxies. Failure to find 
such an offset would strongly suggest that SMBHs grow 
along with disks. 

4.2. Predicted offset from photometric samples 

In order to provide a quantitative prediction for the 
change in the zero-point we use two samples of galaxies 
with detailed photometric d ecompositions. The first is 
the sample of lGadottH ()2009l ) , for which each galaxy is fit 
with three components; a bulge, a disk, and, where nec- 
essary, a bar. We split the sample by whether the galaxy 
is barred or not. In our analysis, for the barred galaxies 
we treat the bar as part of the disk when computing D/B 
values and use Rd from the disk, not the bar. Because 
bars are at the centers of galaxies, our assumption under- 
estimates the fraction of disk+bar mass within the bulge 
effective radius, and therefore also (J 7 ). The second sam- 
ple is the complete and volume-limited catalogue of 86 
low-inclination disk galaxies of all Hubble types observed 
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Fig. 11. — Contours of predicted T. in the B/D-R c / R d plane assuming the [Gultekin ct al. (2009) value of /3 = 4.24 (solid lines). The 
dashed lines correspond to our simulation grid, with horizontal lines at fixed D/B and roughly vertical lines at constant R^ /R e ,Q. A galax y 
growing a disk at fixed R^ evolves from top to bottom parallel to the dashed lines. Left: The (blue) circles are data from Gadotti (2009), 
with open and filled circles corresponding to barred and unbarred galaxies , respect i vely. See text for details of how the parameters for 
the barred galaxies are computed. Right: The (blue) circles are data from Graham (2003), with open and filled circles corresponding to 
galaxies with n < 2 (pseudo) bulge s and n > 2 (classical) bulges, respectively. The (green) star represents the Milky Way based on the 



model of lBissantz fc Gerharc 



bv Ide Jong fe van der KruitJ (|1994f ). For 75 of these, 
iGrahaml (|2003| ) fitted Sersic bulge+exponential disk de- 
compositions in the -ftT-band, regardless of whether they 
are barred or not. We select classical bulges from this 
sample as those galaxies having n > 2, leaving us with 
15 galaxies. 

Figure [11] plots the distribution of both samples in 
the i?d/i? e -D/B plane and overlays contours of T,. All 
but one of the classical bulges in both samples have 
D/B < 1 0, wh ereas many of the pseudo bulges in the 
IGrahaml (|2003P ) sampl e have D/B > 10. The majority 
of the galaxies in the iGadottil (|2009|) sample cluster in 
the range 0.2 < R e /R d < 0.6 (1.6 < R d /R c < 5). For 
about half of a ll galaxie s T, > 1.4, while a small fraction 
(~ 33% ofthe [Graham] (pOOl sample and - 8% in the 
IGadottil (|2009| ) sample) has T. > 2. A small number of 
galaxies fall outside the simulation grid. In calculating 
T. for these galaxies, we extrapolate Eqn. [5] to outside 
our simulation grid. These are mostly however at large 
D/B, and populated solely by pseudo bul ges, rather than 
classical ones. Many more galaxies in thejGadottil ()2009l ) 
sample have D/B < 1 than in the IGrahaml ()200*3l ) sam- 
ple. The difference cannot be attributed to the different 
photometric decompositions since the sa me difference is 
present also for unbarred galaxies in the IGadottil (|2009f ) 
sample. Therefore together these two samples should 
give some indication of the uncertainty in the photo- 
metric parameter (J 7 ). The distributions of T for both 
samples are shown in Figure [12] and the results listed in 
Table [2] Notwithstanding the differences between the 
samples, we find a nar row range of 1 .096 < (J 7 ) < 1.122. 
Assuming /3 = 4.24 (IGultekin et al.l I2009T) . this corre- 



sponds to 1.48 < (r.) < 1.63 and -0 2 1 < 5 a < -0.17. 
For their full sample IGultekin et al.l (|2009f ) measured 
a = 8.12±0.08; with such a small uncertainty on a, off- 
sets between classical bulges and elliptical galaxies in the 
M. — <7 e relation should be measureable. We estimated 
above a scatter due to modelling uncertainties of order 
30%, but the effect we are looking for here is systematic, 
so it should be detectable if present. 



Sample 



N a 



Unbarred Gadotti (20091 
Barred Gadotti ('2009') 
All Gadotti (2009) 
Classical Graham (2003) 



166 
80 

246 
15 



1.098 ±0.004 
1.096 ±0.005 
1.098 ±0.003 
1.122 ±0.025 



TABLE 2 

Values of (J 7 ) for different photometric samples. N g 

INDICATES THE NUMBER OF GALAXIES IN EACH SAMPLE. THE 
ERROR ON (T) IS PURELY STATISTICAL. 



5. TESTING FOR OFFSETS 

We now test for offsets between the M, — a e relations 
of elliptical galaxies and of classical bulges. We show that 
there is no significant offset between the two populations. 
If we decompress the bulges using Eqn. 13 then a small 
but significant offset occurs, which supports our claim 
that an offset should have been detected if SMBHs had 
not grown in mass along with disks. 



5.1. Gultekin sample 
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Fig. 12. — The distribution of T for the various photometric 
samples. The d otted (green) line shows the classical bulg es in the 
IGrahaml [12001) sample. The remaining lines arc for the Gadotti 
l|2009l ) sample of classical bulges: the dashed (red) line is for un- 
barred galaxies, dot-dashed (blue) line for barred galaxies and solid 
(black) line for all galaxies. For ease of comparison, all distribu- 
tions have been normalized to uni t peak value. The to p border is 
labelled by T. assuming = 4.24 HGultekin et al.ll2009h . 



iGiiltekin et al.l <|2009t ) presented a sample of 49 galaxies 
with M, measurements, to which they fitted the M. — a e 
relation. Many of these gala xies have bulge±disk decom- 
posit ions in the literature (jFisher fc Drorvl 120081 120101 
120111) . The ph otometric decompositi o ns we use here 
were t aken from iFisher fc Drorvl (120081 ) , iFisher fc Drorvl 
(|2010f ) and IFisher fc Drorvl (|2011| ). For the unpublished 
decompositions the Appendix provides a description of 
how they were performed. The left panel of Figure H3] 
presents these photometric decompositions with contours 
o f T, from Eq 



n. [9 overl aid. 

IGiiltekin et al.n|2009[ ) found a = 8.23 for ellipticals 
and a = 8.17 for classical bulges, but their definition 
of classical bulges includes the elliptical galaxies. We 
therefore refit the relation to ellipticals and classical 
bulges separately using the code MPFITEX\F| whic h im- 
plements the MPFIT algorithm (|Markwardtll2009l ). We 
firs t fit the M, — o > relat ion of the elliptical galaxies in 
the IGiiltekin~et""aT1 (pOOfll ) sample: IC 1459, M32, M60, 
M84, M87, NGC 821, NGC 1399 (both measurements), 
NGC 2778, NGC 3377, NGC 3379, NGC 3607, NGC 
3608, NGC 4261, NGC 4291, NGC 4459, NGC 4473, 
NGC 4486A, NGC 4697, NGC 5077, NGC 5576, NGC 
5845, NGC 6251, NGC 7052, A1836 and A3565. We ob- 
tain (3 = 4.06 ±0.40, with zero-point a = 8.21 ±0.07 and 
an intrinsic scatter of 0.30. This measurement is in ex- 
cellent agreement w ith the M, — o e relati on of elliptical 
galaxies obtained bv lGiiltekin et all(|2009D using a differ- 
ent fitting method. This M. — a e relation is shown by the 
solid line in the right panel of Figure [T3l We then selected 
the classical bulges to be those having n > 2, leaving us 
with 16 galaxies: NGC 224, NGC 1023, NGC 2787, NGC 
3031, NGC 3115, NGC 3227, NGC 3245 NGC 3585, NGC 
3998, NGC 4026, NGC 4258, NGC 4342, NGC 4564, 
NGC 4594, NGC 4596 and NGC 7457. For this sam- 

6 http://purl.org/mikc/mpfitexy 



pie we measure, from the photometric decompositions, 
(J 7 ) = 1.12 ± 0.03, comparable to the values predicted 
in Table [2] Based on this value and fixing /3 = 4.06, we 
expect 6a = —0.20 if M. had not changed as the disks 
grew. Fitting the M. — a e relation for classical bulges 
while holding j3 fixed, we obtain a = 8.29 ± 0.09, which 
is plotted as the dashed line in FigureQS] The offset from 
the elliptical relation is only 0.08, within the one sigma 
uncertainty and significantly smaller than expected from 
the photometric decomposition if no SMBH growth had 
occurred. The zero-point predicted by the photometric 
decompositions is ~ 2 sigma away from the one found. 
Thus we find no evidence of a significant offset between 
elliptical galaxies and observed classical bulges. 

In order to demonstrate that bulge compression should 
have produced an offset that is measureable, we also 
fitted the M. — <r e relation of the same bulges decom- 
pressed using Eqn. again fixing /3 = 4.06 and using as 
uncertainties on cr ei o and M, the values for the observed 
bulges. We obtain a = 8.47 ± 0.11. The offset from the 
relation for ellipticals is now +0.26, which is two sigma 
different. This fit is shown in the right panel of Figure 
[TBI as the dot-dashed line. 

5.2. Beifiori sample 

As a further demonstration of the absence of an off- 
set in the M. — a e relation between ellipticals and clas- 
sical bulges we cons i der also the in d epend ent sample of 
IBeifiori et all (|2009D . IBeifiori et all (|2009D obtained up- 
per limits on the masses of SMBHs in over 100 galaxies. 
They showed that their relation is parallel to the usual 
M. — cr e relati on, with j3 = 4 . 12 ± .38. For a number of 
these galaxies, IBeifiori et al.l (|2012l ) provide bulge±disk 
decompositions; the resulting sample has 22 disk galax- 
ies. Of these, 16 galaxies have n > 2 which we select 
as classical bulges: NGC 2911, NGC 2964, NGC 3627, 
NGC 3675, NGC 3992, NGC 4203, NGC 4245, NGC 
4314, NGC 4429, NGC 4450, NGC 4477, NGC 4548, 
NGC 4579, NGC 4698, NGC 5005 and NGC 5252. The 
left panel of Figure Q3I plots the distribution of these 
bulges in the B/D-R c /R^ plane. From their photometric 
decompositions we obtain (J 7 ) = 1.16 ±0.03. 

Because this sample only has upper limits on M., not 
actual measurements, we fit the M. — a e relation keep- 
ing the slope of the relatio n fixed to that obtaine d for the 
elliptical galaxies from the IGiiltekin et all (120091) sample, 
i.e. P = 4.06. Although the IBeifiori et al.l (|2009l) sample 
of upper limits cannot give the absolute zero-point of the 
relation, we are interested in relative offsets, for which it 
is well suited. As an estimate for the error on M, we use 
hal f the difference betw een the two upper limits given 
bv IBeifiori eTail (|2009l ). which are based on assuming 
two different inclinations for the nuclear disk surround- 
ing the SMBH. Using a constant error of 10 3 M Q instead 
yields results that are virtually indistinguishable. We fit 
the zero-points for ellipticals (a = 8.46 ±0.10), observed 
classical bulges (a = 8.57±0.10) and decompressed clas- 
sical bulges (a = 8.95±0.11). These results are shown in 
the right panel of Figure [TH The offset between the el- 
lipticals and the observed classical bulges is 0.11, which 
is again less than one sigma. In comparison, the off- 
set between ellipticals and decompressed classical bulges, 
shown in the right panel of Figurc[l3l is ±0.49 (the pho- 
tometric prediction being 5a = 0.26), which is different 
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Fig. 13. — The sample of Giiltckin et al. (2009) with bulgc+disk fits as described in the text. Left panel: Distribution of the sample 
in the B/D-i? /iJ(j plane. Filled symbols show n > 2 (classical) bulges while open symbols show n < 2 (pscudo) bulges. Contours of T. 
assume = 4.24, as before. Right panel: The M. — cr e relation for the 16 classical bulges. The (black) circles with larger values mark 
the observed cr e while the connected (red) squares show a e fi (i.e., the decompressed values) from the fit of Eqn. [9] The solid line shows 
our re-fit for the M. — cr e relation of elliptical galaxies (a = 8.21, /3 = 4.06) while the dashed and dot-dashed lines show the observed 
(a = 8.29 ± 0.09) and decompressed (a = 8.47 ± 0.11) classical bulges fitted by a relation with = 4.06. The shaded region shows the one 
sigma uncertainty on the M. — cr e relation of ellipticals. 

at more than three sigma. Hence, for this sample the ob- 
served offset and the offset predicted if no SMBH growth 
occurs differ by > 3 sigma. 

Thus both the lGultekin et all (120091) sample, with full 
M. measurements, and the iBeifiori et al.l (|2009T ) sample 
with M. upper limits only, show no evidence for an offset 
between the M. — a e relation of ellipticals and of classical 
bulges, even though the disks should have compressed the TABLE 3 

bulges to a mcasureable extent. We therefore conclude Galaxies for which the photometric data predicts r. > 3 

t i ? ai\/l"DTJ 1 • 1 1, 1 1 U • 1 ASSUMING THE M. — a e RELATION OF ELLIPTICALS 

that SMBHs m classical bulges have been growing along ^ p) = (g ^ Am)) The column labelled IReference > 

With dlSKS. LISTS THE SOURCE FOR THE BULGE+DISK DECOMPOSITION. 



NGC 4594 
NGC 3675 
NGC 438 
NGC 3627 
NGC 3140 



6.6 
1.1 
3.1 
3.0 
3.0 



Fisher fc Drorv (20111 
Beifiori et al. (2012) 
Graham (2003) 
Beifiori et al. (2012) 



Graham 



5.3. Galaxies with T. > 3 

We estimated above that the scatter in a e /<7 e .o is 
< 30%. This implies that galaxies in which T > 1.3, 
(i.e. r. > 2.9 for = 4) should be dominated by com- 
pression. Table |3] lists the five galaxies for which the pho- 
tometric properties imply T. > 3; these are the galaxies 
for which the impact of compression is the largest, and 
are therefore ideally suited to test whether or not disk 
(re-)assembly is associated with growth of the SMBH. 
Of the five galaxies, only one, NGC 4594 (the Sombrero 
galaxy), which happens to have the largest T., has a 
proper M, mass measuremen t. Two of the other galax- 
ies have M, upper limits from lBeifiori et al.l (|2009D . The 
remaining two galaxies have no M, measurements that 
we are aware of. We recommend measurements of M, in 
these galaxies in order to further constrain the ability of 
SMBHs to grow along with disks. 

Figure [15] plots the ratio of the observed M. to that 
predicted for its value of a e by the M. — a e relation of 
ellipticals versus r. p hot = J~@, the value of T, predicted 



by the photometric decompositions. The dashed line 
showing M./M. iPro d = r.~ hot represents the location 
of SMBHs that form on the M. — er e relation and do not 
grow as the disk regrows. The dotted line instead shows 
the case M./M. jPre d = 1, corresponding to SMBHs that 
always stay on the M, — a e relation as the disk regrows. 
Most galaxies are above or near the dotted line, and this 
is especially true at T. ph ot > 3, regardless of whether 
the IGultekin et ail (|2C]09T) or the IBeifiori et ail (|2009f ) 
sample is considered. For the galaxy with the largest 
predicted r. iP hot, NGC 4594, we also p lot the improved 
M. measurement of lJardel et al.l (|2011l ) together with its 
uncertainty. NGC 4594 provides the greatest leverage 
in distinguishing how SMBHs and disks co-evolve; Fig- 
ure [15] shows clearly that its SMBH continued to grow 
while its disk was forming. Galaxies would have followed 
the dashed line in Figure [T5] if the M. — Mbui had been 
the more fundamental scaling relation rather than the 
M. — cr e relation, as assumed here. 
There is a hint that barred galaxies are more frequently 
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Fig. 14. — The sample of galaxies with upper limits on M. from Bcifiori ct al. (2009). Left panel: Distribution of the sample in the 
B/D-Rc/ fid plane. Filled sym bols show n > 2 (classical) bulges while open symbols show n < 2 (pseudo) bulges. Contours of T. assume 
/3 = 4.24 HGiiltekin et al.H2009l ), Right panel: The M. — <r e relation of the 16 classical bulges. All black hole masses are upper limits only. 
The (black) circles with larger values mark the observed <r e while the connected (red) squares show cr e ,o (i.e.; the decompressed values) from 
the fit of Eqn. [9] The vario us lines show our fits o f the M. — cr e relation to different samples with the slope of the relation held fixed to that 
for elliptical galaxies in the Giiltckin et al. (2009) sample (/3 = 4.06). The solid line shows the fit to the ellipticals (a = 8.46 ± 0.10), the 
dashed line fits the observed classical bulges (a = 8.57 ±0.10) and the dot-dashed line the decompressed classical bulges (a = 8.95 ±0.11). 
The shaded region shows the one sigma uncertainty on the M. — a e relation of ellipticals. 

found near or below the dashed line in Figure I15[ al- 
though this is not true of all barred galaxies. However 
the data do not reach r # p hot values large enough to de- 
termine whether there is a real difference between barred 
and unbarred galaxies. 

6. DISCUSSION AND CONCLUSIONS 

Observations find th at the peak of the integrated AGN 
activity is at z ~ 2 (jWolf et al.ll2003| ). The majority 
of br i ght quasars are in elliptical galaxies (iKukula et al.l 
[200ltlDunlop et alj|2003t iKauffmann et al.ll2003[) but in- 
termediate brightness Seyfert AGN, which represent a 
significant fracti on of the total AGN number density 
at z = 1.5 — 3 (lUeda et al J 12003ft. a re preferentially in 
disk g alaxies (jSchawinski et al.l 120111 ). iSchawinski et al.1 
(|201lD estimate that 23 - 40% of SMBH growth in these 
AGN occur during a slow, secular mode of the type en- 
visaged here. For the samples of disk galaxies with clas- 
sical bulges that we explored here we estimate a mean 
M. growth by ~ 50% — 65%. If instead we consider the 
samples with M. measurements (both the direct and up- 
per limits only) we find growth factors ~ 60% — 80% 
from their photometric decompositions. The growth fac- 
tor may be somewhat larger still if disks are steeper than 
exponential at their centers. Nonetheless, our estimated 
growth factor spans a range that is broadly in agreement 
with observational estimates. 

We failed to find a significant difference between the 
M. — cr e relation of ellipticals and of classical bulges. 
With currently available samples this result is statisti- 
cally significant only at about two-three sigma. Besides 
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Fig. 15. — The residuals from ( our fit for) the M. — a e rela- 
tion for ellipticals (a = 8.21 for the Giiltckin ct al. (2009) sample 
and a = 8.46 for the IBeifiori et al.l ( 120091) sample, with = 4.06 
for both) versus T. predicted by the photometric decompositions. 
Barred and unbarred galaxies are shown as open an d filled symbols, 
respec tively. The (black) circles are SMBHs from Giiltckin ct al. 
(120091) while T. is compu ted using the bulge+disk decomposition 
of IFisher fc Drorvl (|2010l ) (see text for deta ils). The (red) star s 
are SMBH upper limits fr om the sample of Bcifiori ct al.l 1120091 ). 
with decompositions from IBeifiori et al.l (120121 ). The (blue) trian- 
gle with error bars shows NGC 4594 (t he Sombrero galaxy ) with 
improved M. measurement taken from ljardel et al.l lj2011l ). The 
dashed line shows M./M. , prc d = T. -1 while the dotted line shows 
M./M., pred = 1. 
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increasing the sample size, the best future prospects for 
improving the significance of this result is if more galax- 
ies with photometrically predicted large values of T. were 
to have their M, measured. We have provided a list of 
5 galaxies (Table [3]) with T, > 3; of these 4 have no di- 
rectly measured M.. Galaxies with such large predicted 
growth factors offer excellent probes of the co-evolution 
of SMBHs and disks. Moreover Eqn. [9] makes it easy to 
trawl through photometric catalogs to search for further 
examples of galaxies with large predicted growth factors. 

Since this paper w as first submitted the re have been 
several updates to the lGultekin et all (|2009l ) sample used 
in this work. We exp lored the impact of the se via 
the sample compiled in IMcConnell k Mai (|2012t ). The 
main changes for elliptical galaxies were updates of some 
SMBH masses and <7 e 's and the addition of SMBH mea- 
surements in several brightest cluster galaxies (BCGs). 
Because BCGs evolve differently, we exclude these new 
galaxies from our sa mple and use the sam e sample as 
listed above under the lGultekin et all (|2009l ) sample, up- 
dating to the new M, and a e values (dropping NGC 
2778 which does not have a significant SMBH detec- 
tion in recent measurements). For this sample we obtain 
(a,0) = (8.37±0.07,4.39±0.42). The IMcConnell k Mai 
(|2012t) sample includes a number of new SMBH measure- 
ments in disk galaxies; at present we cannot determine 
whether any of these galaxies host classical bulges. In 
any case, several of these are low mass galaxies and are 
likely to host pseudo bulges, so w e continue to fit to the 
same classical bulge sample from iGli ltekin e t al.1 ()2009l ) 
as before, now fixing the slope to j3 = 4.39. We obtain an 
intercept a = 8.32 ± 0.09, which is statistically indistin- 
guishable from the value for ellipticals. Instead the value 
for the decompressed bulges isa = 8.51±0.11, one sigma 
different from the value for ellipticals. We conclude that 
the latest measurements continue to show no evidence of 
an offset between ellipticals and classical bulges. 

Assuming our result continues to hold with increased 
sample size, the consequence of our finding is that 
SMBHs grow along with disks. The main parameter reg- 
ulating their growth is then the potential within which 
they reside, which is largely set by the bulg e. This means 
that SMBH growth is self-regulated (e.g. iTreister et al.l 
12011 : SMBHs can grow until their feedback unbinds 
any gas otherwise destined to accrete onto them. This 
picture accounts also for the absence of correlations 
with properites of the dark matter halo or of the disk 
(jKormendv k Benderl l201lt iKormendv et al.ll2011f ). 

Disk mass growth leads to an evolution of a e that is 
non-hierarchical, thereby adding nothing to the mass of 
a classical bulge. Another consequence of the absence 
of an offset in the M, — a e relation of classical bulges 
therefore is that the bulge mass, which does not change 
as D/B increases, is not the main parameter determining 
M. . Thus the M. — Mbui relation cannot be as funda- 
mental as the M, — o~ e relation. One interpretation of 
SMBH scaling relations views them as reflecting only a 
central-limit-theorem non- causal evolution produced by 
repea ted galaxy merging (|Pend 120071 : Uahnke k Macciol 
120111 ). In this picture the main correlation is between 
M, an d Mbui, both of which grow during mergers. iPensJ 
(l200l even predicted that bulge-dominated galaxies will 
have tighter scaling relations than disk-dominated ones. 
The lack of an offset between ellipticals and classical 



bulges is contrary to this scenario: some form of reg- 
ulation between SMBHs and bulges is required. 

6.1. The Milky Way Galaxy 

Whether the Milky Way hosts a classical or pseudo 
bulge remains unclear. While its bulge stars are mostly 
old, metal-rich a nd a-enhanced, favoring f ast formation 
during mergers (iMcWilliam k Rich! U995 iZoccali et al l 
2l M I200al2008t iLecureur et all 120071: iFulbright et al.l 
20071 ) , kinematics and morphology f avor its formation 
via the centr al bar (iFul (fl99^H999h : iShen et al.l (pOlOl ) 
but see also iSaha et al.1 (|2012t )). Assuming it is a clas- 



sical bulge, the green star in the right panel of Figure 
ITT1 represents the Milky Way based on a bulge+disk de- 
compo sition of the density model of lBissantz k Gerhard! 
(pOOa ) (D/B = 8.3, R d /R c = 3.2); thi s implies r. ~ 3.7 
If currently M. = 4.1 x 10 6 M Q (IGhez et al.1 120081: 
iGillessen et "aT1l2009h . the original SMBH would have had 
M. - 1.1 x 10 6 M if it formed on the M. - a e relation. 

6.2. Summary 
Our main results can be summarized as follows: 

1. When a disk forms and grows around a pre- 
existing bulge, it gravitationally compresses the 
bulge, causing its effective velocity dispersion, er e , 
to increase. We have provided a fitting formula, 
Eqn. [9j for the change in a e for given bulge-to-disk 
mass and size ratios. 



2. Using the SPSS dat a of IGadottil (120011 and 
iGadotti k Kauffmannl (|200a T~ we find evidence 
that classical bulges have been compressed as disks 
reformed around them. The photometric samples 
predict that bulges should experience a mean in- 
crease in er e by ~ 10%. While small, the steepness 
of the M. — cr e relation requires SMBHs to grow, 
on average, by ~ 50% and extends to > 200%. 

3. The weak correlations between D/B and ag and be- 
tween Rd/Re and 0% ensure that the main effect of 
bulge compression on the M. — a e relation, if M. 
remains unchanged as the disk rcgrows, is an off- 
set to a smaller zero-point at fixed slope. The pre- 
dicted offset between ellipticals and classical bulges 
is measureable with available samples of M,. 

4. We do not find an offset between the M. — <r e re- 
lations of elliptic als and of classical bu lges in ei- 
ther the samp l e of [ Giiltekin et al.l (|2009j ) or that of 
iBeifiori et al.1 (j2009f ). Using available photometric 
decompositions of the galaxies, we show that an off- 
set should have been found if M. had not changed 
since the bulges formed. Thus SMBHs must have 
grown along with disks. 

5. We estimate that SMBHs had to have grown by 
~ 50% — 80% in order to remain on the M. — a e 
relation. Such significant SMBH growth is in agree- 
ment with recent observations that find that at 
1.5 < z < 3 SMBHs in disk galaxies grow by 
~ 23% - 40%. 

6. We have provided a list of 5 galaxies (Table [3]) 
for which the SMBH is predicted to have needed 
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to grow by a factor greater than three to remain 
on the M. — a e relation. SMBHs with such large 
growth factors provide strong constraints on the 
mechanisms regulating the M, — a e relation and 
we strongly encourage measurement of their black 
hole masses. 
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APPENDIX 

UNPUBLISHED PHOTOMETRIC FITS 

We make use of unpublished decompositions of disk galaxies in the sam ple oflGiiltekin et al.l (l2009lj kindl y provided 
to us by David Fisher. Most of these decompositions have been published (jFisher fc Drorvll2008ll2oToll2011| ). However 
a few remain unpublished and we provide here a description of the analysis method by which David Fisher derived 
these decompositions. 

The decompositions use archival HST and groun d-based data. When p ossible, near infrared data are used as they 
are less sensitive to the obscuring effects of dust. iFisher &: Drorvl (|2008| ) show that for relative quantities, such as 
B/T, there is little difference from U-band to iJ-band; data are therefore restricted to be U-band or redder. For 
each galaxy the surface brightness profile is determined through ellipse fitting of both HST and ground-based data, 
thereby simultaneously constraining both the small scale structure at the center of the galaxy and the shape of the 
outer disk profile. Interfering objects, such as foreground stars and background galaxies, are masked via automatic 
source identification methods and manually removed. For ground-based images the sky is removed by subtracting a 
surface, fitted to regions of images that do not contain galaxy light. The radial sizes of the ellipses are optimized to 
maintain a roughly constant signal-to-noise ratio across the profile, and zero-point shifts of the ground-based image 
to match the HST data ensured continuity. The bulge+disk decompositions are then determined by fitting a Scrsic 
bulge plus outer exponential disk to the major axis surface brightness profile. 
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